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A synthesis of Polyakov loop models of the deconfinement transition and quasiparticle models of gluon plasma 
thermodynamics leads to a class of models in which gluons move in a non-trivial Polyakov loop background. These 
models are successful candidates for explaining both critical behavior and the equation of state for the SU(3) 
gauge theory at finite temperature. Two sensitive lattice probes which can discriminate between models are the 
Polyakov loop P and e — 3p. 



The body of results on the critical behavior and 
thermodynamics of QCD and related theories at 
finite temperature are among the most important 
that lattice gauge theory simulations have pro- 
vided, giving important guidance in both theory 
and experiment. Theoretical work on the ther- 
modynamics of the quark-gluon plasma at zero 
chemical potential is often benchmarked against 
results from lattice simulations. Here we explore 
a class of models which synthesize the attrac- 
tive features of quasiparticle models of the quark- 
gluon plasma and Polyakov loop mod- 
els of the deconfinement phase transition 00. 
Quasiparticle models, which are partially moti- 
vated by perturbation theory, have been success- 
ful in fitting lattice results for thermodynamic 
quantities such as the pressure over a large range 
of temperatures, but typically at the cost of in- 
troducing a number of extra assumptions and 
parameters. Polyakov loop models extend the 
equilibrium free energy density / to a function 
of the Polyakov loop P, which is the order pa- 
rameter for the deconfining phase transition in 
pure gauge theories. The spontaneous breaking 
of center symmetry in pure gauge theories, which 



underlies our understanding of the deconfinement 
transition, occurs in these models because the free 
energy develops symmetry-breaking global min- 
ima above the deconfinement temperature T c . 

The general class of models of deconfinement 
we study have free energies of the form 



f(T,P) = V(T,P)-p g (T,P,M(T)) 
+B(T,P,M(T)). 



(1) 
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The first term V(T, P) is phenomenological, and 
favors the confined phase at low temperature. 
The second term is the negative of the quasipar- 
ticle pressure p g , and represents the contribution 
of gluonic quasiparticles moving in the presence 
of a background Polyakov loop. The quasiparticle 
pressure is given by 
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■Tr A In [1 - e- 0u ">P] (2) 

where TrA denotes the trace in the adjoint 
representation. The pressure depends on 
the quasiparticle mass through the quasiparti- 
cle energy, assumed to have the simple form 

u>k = \J k 2 + M (T) 2 ; the quasiparticle mass is 
temperature-dependent. The third term is the 
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so-called bag term. It arises because the zero- 
quasiparticle state has a non-zero, temperature- 
dependent energy B. If we assume that B de- 
pends on T only through the quasiparticle mass 
M, then thermodynamic consistency leads to an 
expression for B which may be written as 



B(T,P,M(T)) 



T dM 2 f d 3 k 1 
To T ~dr~ J (2tt) 3 k 2 

■Tr A In [l - e-^^p] . (3) 



After minimization with respect to P, the equi- 
librium pressure is given by the negative of the 
minimum of the free energy density. 

In previous work 0, we have developed two 
models for the SU(N) deconfinement transition 
which fit within this general class. In both mod- 
els, we took M = and hence B = 0. Both of 
these models gave a good representation of the 
critical behavior for all N, and reasonable ther- 
modynamic behavior. Here we extend the second 
of these two models, called model B, to model 
B+M which includes quasiparticle mass effects. 
The potential V associated with model B has the 
form 



V(T,P) 
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where p(P) is Haar measure on the gauge group. 
This form for the potential can be understood as 
enforcing color neutrality at scales R and larger. 
We again use the parameter R to set T c to the 
result of lattice simulations. The constant Vq is 
used to match the pressure at one point on the 
curve, again chosen to be T c . 

We show in figure ^ the dimensionless pres- 
sure p/T 4 as a function of T/T c for models B and 
B+M, along with simulation data from Note 
that once R and vo are fixed by determining T c 
and p(T c ), there are no more free parameters in 
model B. The pressure from Model B overshoots 
the lattice data at higher temperatures, and the 
quasiparticle mass corrects this. The quasipar- 
ticle mass term has the form suggested by per- 
turbation theory, M(T) = g{T)T/y/2. We take 
the temperature-dependent coupling constant to 
be given by g 2 {T) = 8vr 2 /H In (T/A) where A is 
an adjustable parameter. We have obtained the 



best fit to the lattice pressure data by replacing 
T/A in g(T) with a crossover function of the form 

[(T/A)" + (T cr /A) n ] 1/n where we have taken n = 
10. Both Polyakov loop effects and a rapidly ris- 
ing quasiparticle mass suppress the temperature 
as one moves to lower temperature, and both are 
not needed. It is physically plausible that g does 
not run with T all the way down to zero temper- 
ature, and the crossover function provides that 
effect. The critical temperature is taken to be 
T c = 0.272 GeV, although an exact number is not 
needed. The parameter R is determined from T c , 
and is given by R- 1 = 0.61T C = 0.17 GeV which 
may be compared with the result for model J5, 
R- 1 = 73T C = 0.2 GeV. The value of A obtained 
was 0.3T C = .081 GeV, with a crossover tempera- 
ture T„ = 2AT r . 
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Figure 1. p/T 4 versus T/T c . 

Figure El compares the dimensionless inter- 
action measure A = (e — 3p) / T 4 for model 
B+M with both lattice results and a quasiparti- 
cle model without Polyakov loop effects (L — 1). 
Model B+M shows a peak in A just above T c , 
consistent with the lattice data, and correctly re- 
produces the latent heat of the first-order SU (3) 
deconfinement transition. Of course, the model 
without Polyakov loops effects is incapable of 
reproducing the deconfinement transition. The 
first-order character of the SU(3) deconfinement 
transition is also clear in Figure OH which shows 
the Polyakov loop as a function of T/T c The shape 
of the curve is similar to that found in a recent 
lattice determination of the Polyakov loop [7]- It 
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is also interesting to compare the quasiparticle 
masses for model B+M with the L = 1 model, 
as shown in figure 01 The most notable difference 
is the rapid rise of the quasiparticle mass in the 
L = 1 model near T r . 
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Figure 2. (e - 3p)/T 4 versus versus T/T c 
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Figure 3. L = Tr F (P) /N c versus T/T c . 

It should be apparent that many models will 
be able to fit lattice thermodynamic results well. 
The ability to describe both the thermodynamic 
results and the deconfinement transition, while 
not unique to the model described here, seems 
very desirable. The introduction of an extra term 
V in the free energy is phenomenological, but 
based on our understanding of the deconfinement 
transition. The onus of describing the tempera- 
ture region just above T c is placed where it likely 
belongs, on the confinement mechanism. It is pos- 
sible to extend this model to include quarks. As 
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Figure 4. M(T)/T versus versus T/T c . 



with confinement, our inability to describe chiral 
symmetry breaking in a fundamental way forces 
us to use phenomenological models, in this case of 
Nambu-Jona Lasinio type. Nevertheless, a useful 
synthesis of our current understanding of decon- 
finement and chiral symmetry restoration is likely 
possible. 
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